Goos-Hanchen-like shifts for Dirac fermions in monolayer graphene barrier 
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We investigate the Goos-Hanchen-like shifts for Dirac fermions in transmission through a mono- 
layer graphene barrier. The lateral shifts, as the functions of the barrier's width and the incidence 
angle, can be negative and positive in Klein tunneling and classical motion, respectively. Due to their 
relations to the transmission gap, the lateral shifts can be enhanced by the transmission resonances 
when the incidence angle is less than the critical angle for total reflection, while their magnitudes 
become only the order of Fermi wavelength when the incidence angle is larger than the critical angle. 
These tunable beam shifts can also be modulated by the height of potential barrier and the induced 
gap, which gives rise to the applications in graphene-based devices. 

PACS numbers: 73. 23. Ad, 42.25.Gy, 72.90.+y, 73.50.-h 



I. INTRODUCTION 

Monolayer graphene has attracted much attention-^ 2 - 
since the graphitic sheet of one-atom thickness has 
been experimentally realized by A. K. Geim et al. in 
2004 3 . The valence electron dynamics in such a truly 
two-dimensional (2D) material is governed by a mass- 
less Dirac equation. Thus graphene has many unique 
electronic and transport properties^ 2 -, including Klein 
tunneling 4 . Recently, further investigations show that 
a method to produce a finite bandgap in graphene sheets 
by epitaxially on proper substrate^ has been proposed, 
and the Dirac fermions in gapped graphene are described 
by 2D massive Dirac equations. The induced gap at the 
Dirac point is significant to control the transport of the 
carriers and integrating graphene into the semiconductor 
technology. 

Motivated by these progress, the transport of massive 
or massless Dirac fermions in graphene opens a way to 
design various graphene-based electron devices in term 
of the electron optics behaviors^—, such as focusing 6 , 
collimation 7 , Bragg reflection^, and Goos-Hanchen effect 
(GH) 10 ' 11 . In particular, Zhao and Yelin 1 ^ have shown 
that, based on the electronic counterpart to the trapped 
rainbow effect in optics^ 2 -, the interplay GH effect and 
negative refraction 6 - in graphene leads to the coherent 
graphene devices, such as movable mirrors, buffers and 
memories. Beenakker et alr^ have further found that the 
GH effect at a n-p interface in graphene doubles the de- 
generacy of the lowest propagating mode, which can be 
observed as a stepwise increase by 8e 2 /h of the conduc- 
tance with increasing channel width. 

In this paper, we will investigate the negative and 
positive lateral shifts for Dirac fermions in transmission 
through a 2D monolayer graphene barrier, based on the 
tunable transmission gap 1 -. Generally, the magnitude of 
GH shift for total reflection in graphene is in the order of 
Fermi wavelength. However, the lateral shifts discussed 
here are similar to but different from the conventional 
GH shift, because they do result from the beam width 
of Dirac fermions, and can be enhanced by the transmis- 



FIG. 1. (a) Schematic diagram for two-dimensional mono- 
layer graphene barrier, (b) Negative and positive lateral shifts 
of Dirac fermions in transmission. 



sion resonances. As a matter of fact, the lateral shifts 
in transmission have nothing to do with the evanescent 
wave, thus we term them as Goos-Hanchen-like (GHL) 
shift, which can be considered as an electronic analog of 
the lateral shifts in the optics^ 4 , and atom optics 1 ^. More 
interestingly, not only large positive but also large nega- 
tive GHL shifts can occur in transmitted beam through 
the graphene barrier. The negative shift behaves like the 
phenomena of negative refraction in graphene 6 -. In addi- 
tion, we will also discuss the effect of the induced gap on 
the GHL shifts in the gapped graphene barrier. All these 
tunable beam shifts in the monolayer graphene barrier 
can be applied in the design of graphene-based devices, 
such as electron wave switch, wave vector or energy filters 
and splitter. 



II. THEORETICAL MODEL 

Consider the massless Dirac fermions with Fermi en- 
ergy E at angle c£>q with respective to the x axis incident 
from zero-gap graphene upon a 2D gapped graphene bar- 
rier of hight Vq in which the electron acquires a finite 
mass of A/u 2 , as shown in Fig. [1] where the tunable 
potential barrier is formed by a bipolar junction (jhn-p) 
within a single-layer graphene sheet with top gate volt- 
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age V^, Vq and d are the height and width of poten- 
tial barrier, respectively. Since graphene is a 2D zero- 
gap semiconductor with the linear dispersion relation, 
E = HkfVf, the massless electrons are formally described 
by the Dirac-like hamiltonian^, Hq = —ihvfaV, where 
Vf « IC^to-s^ 1 is the Fermi velocity, kf is the Fermi wave 
vector (A = 2ix jkf is Fermi wavelength), and a = (cr x , a y ) 
are the Pauli matrices. The wave functions of the plane 
wave components for the incident and reflected regions 
are assumed to be 



*z=(^)e^^>+r(_ ae 1 . 
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so the corresponding wave function in the transmitted 
region can be expressed by 



= t 
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where s = sgn(E), k x = kf cos <p and k y — fc^sin0 
are the perpendicular and parallel wave vector compo- 
nents outside the barrier, <fr is the incidence angle of the 
plane wave component under consideration. For the gen- 
eral case of the gapped graphene barrier, the Hamilto- 
nian of massive Dirac fermions can be written down as 
Hi = —ihvfaV + Acr z , where A is equal to the half of 
the induced gap in graphene spectrum and positive (neg- 
ative) sign corresponds to the K (K') point, thus the 
wave functions in the barrier region have the following 
form: 
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where s' = sgn(E - V Q ), k' f = a/(V - Ef - A 2 /Hv f , 
Qx = (kf—ky) 1 ' 2 , 9 = &TCta,n(k y / q x ) , a and (3 are defined 

by 
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Accordingly, the critical angle <j) c for total reflection can 
be defined by 



arcsm 



(V - Ef - A 2 
E 2 
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so that when <j) > (f> Cl the wave function in the propagat- 
ing case becomes evanescent wave by replacing q x with 
in, where k = (k 2 — k 2 ) 1 / 2 . According to the bound- 
ary conditions, the transmission coefficient t = e lv j f is 
determined by 

t _ 1 

cos(q x d) — i(ss'x sec <f> sec 9 + tan <p tan 9) shx(q x d) ' 

(5) 




FIG. 2. (Color online) Dependence of GHL shifts in the 
propagating case on the barrier's width d, where <j>o = 25°, 
Vo = 120meV , d is re-scaled to k x ^d, E = 80meV, A = 
20meV (solid line), E = 220meV ', A = 20meV (dashed line), 
E = 80meV, A = OmeV (dotted line), and E = 220meV, 
A = OmeV (dot-dashed line). 



where \ = J A 2 + kf / k'j and the phase shift if is ob- 
tained by 

[ sin^ sin </> + ss'x , , „1 , c . 
ip = arctan - tan (q x d) . (6) 
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From the above expression, it is clear that in the limit 
A — > 0, we get \ — 1 and thus one can obtain the same 
expression for electronic transmission probability T = 
l// 2 corresponding to the massless Dirac fermions^. 



For a well-collimated beam with the central angle <fio 
of incidence, the GHL shift can be defined, according to 
the stationary phase metho d 10 ' 17 , as 
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where the subscript in this paper denotes the values 
taken at <fi — 4>q. The most intriguing behavior in prop- 
agating case is found for Klein tunneling, E < Vq, where 
the GHL shifts can be negative, and also be enhanced by 
the transmission resonances, whereas the lateral shifts for 
classical motion, E > Vq, are always large and positive. 
Examples of their dependence on the barrier's width at 
different induced gaps A are plotted in Fig. [21 where 
E = 80meV, V = l20meV, (j) = 25° is less than the 
critical angle defined by Eq. (j4]), solid and dotted lines 
correspond to Klein tunneling E = 80meV < Vo, and 
dashed and dot-dashed ones correspond to classical mo- 
tion E = 220meV > Vq. On the contrary, when the 
incident angle (j)Q is larger than the critical angle <j) c , the 
lateral shifts become in the order of Fermi wavelength due 
to the evanescent wave, which is similar to those in total 
reflection at a single graphene interfac o 10 ' 11 . Instead of 
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FIG. 3. (Color online) Dependence of GHL shifts in the 
evanescent case on the barrier's width d, where (a) <j>o — 35° 
and (b) <j>o = 75°, d is re-scaled to Kod, E — 80meV , A = 
20meV (solid line), E = 220meV, A = 20meV (dashed line), 
E = 80meV, A = OmeV (dotted line), and E = 220meV, 
A = OmeV (dot-dashed line). 



the enhancement by the transmission resonances shown 
in Fig. 03 Fig. |3] illustrates that the GHL shifts for Klein 
tunneling and classical motion saturate respectively to 
negative and positive constants with increasing the bar- 
rier's width in the evanescent case, where (a) cf> — 35° 
and (b) 75°, d is re-scaled to /tod, the other parameters 
are the same as in Fig. [5] 



The GHL shifts obtained above can be positive 
as well as negative, depending on the influence of 
sm(2q x0 d) / (2q x0 d) . Since the inequality 
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the lateral shifts can be positive only when 
sin(2q x0 d) / (2q x0 d) — > 1 for a very thin barrier, 
namely d — > 0. However, the lateral shifts turn negative 
with increasing d. It is interesting that the negative 
lateral shifts can be enhanced by the transmission 
resonances. Since the transmission coefficient is an 
oscillating function of tunneling parameters and can 
be exhibit any value from to 1—, there are resonance 
conditions q x od = Ntt, N = 0, ±1, ±2, ... at which the 
barrier is transparent, T = 1. At resonances, the lateral 
shifts in this case reach s\ qx0 d=NTT = — (&o/<7xo)^tan0O) 
which correspond to the maximum absolute values. At 
anti-resonance, q x od — (N + 1/2)tt, the lateral shifts 
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exotic behaviors of negative and positive GHL shifts are 
analogous to those of lateral shifts for the transmitted 
light beam though left-handed metamaterial slab^*, 
based on the link between Klein paradox and negative 
refraction^. 

On the contrary, when the incidence angle is larger 
than the critical angle <f/ c , <f>o > <p' c , the lateral shift be- 
comes 
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III. DISCUSSIONS 



In the limit of opaque barrier, n d — > oo, the lateral shift 
trends to a constant as follows, 



In this section, we will shed light on the properties of 
GHL shifts in details. For simplicity, we will find the 
following analytical solutions in the limit A = and 
discuss the negative and positive lateral shifts in cases 
of Klein tunneling and classical motion, respectively. 

Case 1: Klein tunneling (ss' = —1). In this case, the 
critical angle ^ becomes 



Vo 
E 



(8) 



when the condition E < Vq < 2E is satisfied. When 
the incidence angle <fio is less then the critical angle <f>' c 
obtained above, 4>q < <p' c , the lateral shift is given by 
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where fco = (^/^/ 
T is given by Eq. © 



kyo) 1 / 2 , and transmission probability 
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which is proportional to 1/k q , and implies that the 
GHL shift in the evanescent case is in the same or- 
der of electron wavelength as the GH effect in a single 
graphene interfac o 10 ' 11 . More interestingly, the saturated 
GHL shift is negative when the incidence angle satisfies 



< 4>o < 4>* , where the critical angle is defined by 
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(12) 



But the GHL shift in this case will becomes positive 
when 4>o > <fi* . The sign change of GHL shifts described 
by Fig. [3] (b) appears at the incidence angle <fio = <fr* , 
which is similar to the result of the quantum GH ef- 
fect in graphene, taking the the pseudospin degree into 
account^. 

Case 2: classical motion (ss' = 1). In this case, the 
critical angle is 



/o 2 



cos 2 \q x od) 



k 4 

-2—2- sin (q x0 d) 

K-xOlxO 



(9) 



in . I , Vo 

b„ = arcsm | 1 — ^- 



(13) 



4 



20 



10 



I- 0.5 



' 1 / 
/*'"*' ^-J" " 






(b) 


P.': 1 ' •". ■ -.1 

>: • V 1 
i -, , i 
1 : | i / 
1 : ; i ' 

1 : 1 

« • i I 

Vs J.J J 





50 100 150 200 250 

E (meV) 



dotted lines correspond to </>o = 20°, 15°, and 10°, re- 
spectively. It is shown that the GHL shifts is closely 
related to the transmission gap AE — 2hk y Vf 13 . Fig. 
2] indicates that the lateral shifts change the sign near 
the Dirac point E = Vq, and can also be enhanced by 
the transmission resonances near the boundaries of en- 
ergy gap. In addition, the incidence angle has also great 
impact on the GHL shifts. The absolute values of the 
lateral shifts increase with increasing the incidence an- 
gles, and the positions of the maximum (absolute) values 
for the positive (negative) beam shifts can also be tuned 
because of the resonance conditions, as shown in Fig. 



FIG. 4. (Color online) GHL shifts (a) and transmission gap 
(b) as the function of the incident energy E, where d = 80nm, 
Vo — 120meV , and A = OmeV . Solid, dashed, and dotted 
lines correspond to cj>o = 20°, 15°, and 10°, respectively. 



When the incidence angle is less than the critical angle for 
total reflection, 4>q < <j>", the lateral shift can be written 
as 
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where k' Q = (fc/fcy — ky^) 1 / 2 , and transmission probability 
is 
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Similarly, the lateral shifts for classical motion also de- 
pend periodically on the barrier's width, thus can be 
enhanced by the transmission resonances. The lateral 
shifts at resonances reach s 



q^od—Nrt 



while at anti-resonance they become s\ qsc0 d=(N+i/2)-n = 
(k xQ /k^)dtari(j)Q. However, these GHL shifts in clas- 
sical motion are always positive as those in the two- 
dimensional semiconductor barrier 17 . When tfio > <fi", 
the GHL shift in the evanescent case becomes 
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Then the lateral shift in the limit, Kgd 
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FIG. 5. (Color online) Dependence of GHL shifts on height 
Vo of potential barrier (a) and induced gap A (b), where d = 
80nm, (a) A = OmeV, E = 150meF, cj> = 20° (solid line), 
E = 150meF, O = 10° (dashed line), E = 100meF, cj> = 
20° (dotted line), E = lOOmeF, 0o = 10° (dot-dashed line); 
(b) Vo = UOrneV, E = 220meV", O = 20° (solid line), E = 
220meF, (j> = 10° (dashed line), E = 80meF, cj} = 20° 
(dotted line), E = 80meF, O = 10° (dot-dashed line). 



which is always positive constant. 



Based on the properties in two cases of Klein tunneling 
and classical motion, the GHL shifts (a) and correspond- 
ing transmission probabilities (b) as the function of inci- 
dence energy E are shown in Fig. 01where d = 80nm, 
Vo = l20meV, and A = OmeV. Solid, dashed, and 



Finally, we will turn to discuss the modulations of GHL 
shifts by the height Vo of potential barrier and the in- 
duced gap A. It is shown in Fig. [5] that the GHL shifts 
can be controlled by changing the height Vb of potential 
barrier, which can be easily implemented by applying a 
local top gate voltage V g to graphenei^. Since the lat- 
eral shifts are in the forward and backward directions in 
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the cases of E > Vq and E < Vq, respectively, it is sug- 
gested that the incidence energy can be selected by these 
tunable beam shifts. Thus, this phenomenon does result 
in an alternative way to realize the graphene-based elec- 
tronic devices, for example, energy splitter and energy 
hltcr. Moreover, Fig. [5] (b) further investigate how the 
GHL shifts are affected by a gap opening at the Dirac 
points. Comparisons of Figs. [Hand Figs. [3] further show 
that the energy gap will increase (decrease) the absolute 
values of the shifts in the propagating (evanescent) case. 
And the influence of gap in the propagating case is more 
pronounced than that in the evanescent one. The method 
to generate the energy gap in graphene is through an in- 
version symmetry breaking of the sublattice due to the 
fact the densities of the particles associated with the on- 
site energy for A and B sublattice are different^. There- 
fore, the periodical dependence of GHL shifts on the gap 
provides an efficient way to modulate the lateral shifts in 
a fixed graphene barrier. 

IV. CONCLUSION 

In summary, we have investigated the GHL shifts 
for Dirac fermions in transmission through a monolayer 
graphene barrier. The lateral shifts, as the functions of 
the barrier's width and the incidence angle, can be nega- 
tive and positive in Klein tunneling and classical motion, 
respectively. Since the lateral shifts have a close rela- 



tion with the transmission probability, the lateral shifts 
can be enhanced by the transmission resonances when 
the incidence angle is less than the critical angle for to- 
tal reflection, while their magnitudes arc only the order 
of Fermi wavelength when the incidence angle is larger 
than the critical angle. Compared with the smallness 
of conventional GH shift in graphene, the large nega- 
tive and positive GHL shifts, which can also be modu- 
lated by the height of potential barrier and the induced 
gap, will have potential applications in various graphene- 
based electronic devices. We further hope that these sim- 
ilar phenomena in magnetic graphene barrier may lead to 
the graphene-based spintronic devices on spin filter and 
spin beam splitte r 17 ! 19 . 
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